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This paper describes a new magnetic trap for ultra-cold neutrons (UCNs) 
made from a 1.2 m long Halbach-octupole array of permanent magnets 
with an inner bore radius of 47 mm combined with an assembly of su¬ 
perconducting end coils and bias field solenoid. The use of the trap in 
a vertical, magneto-gravitational and a horizontal setup are compared in 
terms of the effective volume and ability to control key systematic effects 
that need to be addressed in high precision neutron lifetime measurements. 


1 Introduction 

The free neutron undergoes /3-decay via n —1 p + e~ +D e . Precise measurements 
of the mean lifetime r n are used for obtaining the universal weak coupling con¬ 
stants of the nucleon from which one derives important semi-leptonic weak cross 
sections. They are also needed for searches of beyond Standard Model physics, 
and for calculations of primordial helium abundance in Big Bang Nucleosyn¬ 
thesis. These motivations are described in more details elsewhere in these pro¬ 
ceedings, as well as in various review papers on the neutron particle physics 
field jH[2[ or on r n specifically pH®- 

Ultra-cold neutrons (UCNs) are free neutrons with kinetic energies less than 
the neutron optical potential of well-chosen materials so that they can be con¬ 
fined in a material “bottle” via total internal reflections. For instance, beryllium, 
a commonly-used material for UCN reflection, has Vee = 252 neV, correspond¬ 
ing to a velocity of ss 7ms -1 . UCNs stored in bottlci0for measuring r n have 
been used in the most precise experiments to dated Elans mi Eg. 

Due to the size of the neutron’s magnetic moment /x n of 60neVT _1 , the 
Stern-Gerlach force can be used to confine UCNs also. This technique [131 HU US] 
offers a method for confining UCNs free from energy-dependent wall losses that 
require corrections. While the traditional “counting the survivors” scheme can 

2 Currently at North Carolina State University, 2401 Stinson Dr., Riddick 421, Raleigh, NC 
27695, USA. 

x The exceptions being in-beam measurements [5] EH- 
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be used for determining r n , extracting and detecting the charged decay products 
is also possible with magnetic traps. A strong motivation for using magnetic 
trapping is to improve upon material bottle experiments with a better control 
and reduction of systematic corrections [Tip, 171 T8, T9], 

In magnetic bottles losses of UCNs due to depolarization at the weak field 
regions has been shown theoretically to be suppressible if a sufficiently strong 
bias field is used [231 GSl- However, the effects caused by phase-space evolution of 
the UCN gas and the issue of how to effectively remove above-threshold UCNs 
must be carefully addressed by high-precision measurements. There is also the 
possibility of a gradual warming of the UCN spectrum, caused by magnetic 
noise or mechanical vibrations, which has largely been unexplored. With these 
issues in mind, two concepts using an Halbach-octupole array combined with a 
superconducting coil assembly are discussed in this paper. 

Previously, the idea of using a UCN production volume integrated inside a 
horizontal, sliding magnetic bottle was presented (Hi. Its goal was to extract 
produced UCNs to vacuum in order to avoid losses due to interactions with the 
superfluid helium converter and also to avoid the dilution of the high density 
of UCNs offered by super-thermal production [32] • In this scheme, due to the 
geometry of available cold neutron beams, the long axis of the trap is required 
to be horizontal. However, after further exploration of the idea, and combined 
with the recent success of a high-density superfluid helium UCN source using a 
vertical, window-less extraction system where transmission losses through win¬ 
dows are eliminated[23]. we decided to employ a more traditional scheme. 

2 Design of magnetic fields 

The central component of the discussed setups is the 1.2 m long 32-piece Halbach- 
typeEH octupole array for radial UCN confinement (shown in Fig. [I]). The 
octupole has an inner bore radius R = 47 mm and a nominal B = 1.3 T at its 
surface (at room temperature). While the field near the center is that of an 
ideal octupole B(p) oc p 2 3 , deviations appear near R due to the discrete number 
of magnets. The flux density B from 2D finite element calculations with FEMM 
[25| at different azimuthal angles <f> and distances from R is shown in Fig. [2j 
The weakest B for a fixed radial position p occurs at the off-pole piece^J 

These calculations also show that the magnetic material at the off-pole 
pieces’ inner corners are the most susceptible to demagnetization. To model 
how this might affect B in the bore, these corners were rounded out with a 
1.5mm radius (and replaced with an air gap). This causes sharp dips right at 
the corners and a slight enhancement at the center of the piece (also in Fig. [2]). 
These features become insignificant at distances 2 mm from R. Conclusions 
are the same for different corner geometries (quarter-circular and chamfer cut). 
Besides from demagnetization, if chips in the brittle NdFeB material exist, then 

2 The pieces are called off-pole when they have their magnetization vector parallel or anti¬ 

parallel to the <j> unit vector at the centre of that piece. 
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Figure 1: The octupole magnet array assembly: (a) rods used to fix the rotation 
of the 12 modules, (b) split aluminum shell for holding the modules together, 
(c) stainless steel end plates fixed to the aluminum shell, and (d) brass pressing 
screws. The ribbed structure of the aluminum shell provides good thermal 
contact for cooling the magnets to ~ 120 K. An individual module is shown in 
Ref. [2T]. 



Figure 2: 2D finite element calculations of the Halbach-octupole. B is plotted 
for cj) spanning a quadrant in the cry-plane and for different distances from R 
(indicated by the boxed numbers). The (j> positions of the centers of the pole 
and off-pole pieces are indicated. The fainter lines for each wall distance are 
from modeling the demagnetized magnets with rounded corners. 
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Figure 3: The field at the ends from solely the octupole array: B at R becomes 
weaker and strong B z components start to emerge. A quadrant slice of B z at 
z = 575 mm, the central position of the end coil placed closest to the end of the 
array (see Fig. [4]). 


similar dips could exist. It is therefore wise to keep UCNs from exploring too 
close to these regions. 

To understand the field at the ends of the array 3D finite element calculations 
were performed with R ADI A [25]. This revealed that B(R) drops to < 0.8 T at 
the end of the array and there exist axial components B z « 0.8 T that are 
strongest near the surface of the pole pieces. These features are shown in Fig. [3] 
and have been confirmed with Hall probe measurements [27]. 

A superconducting coil assembly (depicted in Fig. [4]) consisting of a small end 
coil, a bias field solenoid, and a large end coil is used for axial UCN confinement 
and removal of the low-field region. The fields that can be produced by these 
coils when they are run at their maximum current of 300 A are 1.7 T, 1.2 T, and 
510 respectively. The 30 cm inner diameter of the coils produces only small 
B p components at R, and hence there is only a small cancellation with B p from 
the octupole array. The end coils are slotted inside the octupole array so that 
the region of strong B z cancellation is situated away from the trapped UCNs. 
Two magnetic field configurations from the coils that will be discussed in the 

3 This allows focussing decay protons onto a small detector. This might be done in later 
experiments but will not discussed in this paper. 
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Figure 4: Top: a to-scale drawing of the octupole array and the surrounding 
assembly of superconducting coils. Bottom: the axial field from the coils for the 
two trapping configurations. The vertical solid lines show the end of the array 
and dash-dotted lines the center of the end coils. 


next section are also shown in Fig. [4] The bias field is chosen to be > 0.1 T for 
preventing depolarization in the low field region. 

Combining the 3D calculations with all field sources, it was shown that 
demagnetization of the permanent magnets can be avoided by cooling them to 
~ 120 K, which leads to a strong increase in the coercivity. This cooling also 
increases B at R by ~10%, which is not included in the discussion to keep the 
estimates of the trap depth conservative. 

3 Trapping potential and effective volume 

The total energy of a neutron E n at a point r = (x, y , z) in space is given by: 

E n = E kin (r) + E pot (r) = E kin (r) + V glav (z) + V ma g (r) 

= E kin (f ) + mgz ± y n B(r ) , (1) 

where E k i n is the kinetic energy, m is the neutron mass, g is the gravitational 
acceleration, and z is the height. A constant offset is neglected to place Spot at 
the potential minimum of the trap. In the last term, the plus sign refers to the 
low field seeking spin-state. Maps of S pot for a horizontal configuration, where 


5 




















z - horizontal [cm] 


Figure 5: A map of E pot for the horizontal configuration with the axial fields 
shown in Fig. [4] and for the <j> = 90° slice aligned vertically. The thick line 
shows the contour for Strap = 63neV and the dotted contour lines are placed 
at increments of 5 neV. 


both end coils are required for axial confinement, and for a vertical configuration, 
where only the bottom end coil is required due to the gravitational potential 
mg = 102neVm _1 ( magneto-gravitational trap), are shown in Figs. [ 5 ] and [6j 
In order to populate the trap with UCNs, current in the small end coil will be 
lowered temporarily to reduce E pot . 

The depth of a pure magnetic trap S trap is defined by the maximum E n 
a UCN can have without the possibility of it exiting the trap at the ends or 
making contact with material walls. Since an inner tube will be placed inside 
the magnet bore, a reduction in the radius of the trap of 1.5 mm is taken into 
account for determining S trap . 

A UCN clearly cannot explore regions in a trap where E n < S pot (r) so 
that the accessible volume is energy-dependent. Thus, it is useful to define the 
concept of the effective volume[25] V e g{E n ), which can be expressed as: 


V e ff(E n ) = Re 




dV 


( 2 ) 


where f v is over volume of the trap and dV is the volume element at r. This is 

shown for the two configurations in Fig. [71 The total number of UCNs stored 

E ^ 

in the bottle is then given by: f Q trap n(E n ) V e g(E n ) dE n , where n(E n ) is the 
energy-dependent spectral UCN density. If this is calculated for the typical 
Maxwellian spectrum, where n(E n ) oc \/E n , then the number that can be stored 
in the horizontal configuration is ~2.5 times greater than in the vertical. 
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Figure 6: A map of E pot for the vertical configuration with the axial fields shown 
in Fig. [4] for tf> = 0° (maximum B p reinforcement) and <f> = 45° (maximum B p 
cancellation) slices. The thick solid line is for E tlap = 48neV and the dotted 
contour lines are placed at increments of 5 neV. 
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Figure 7: The UCN energy-dependent effective volumes for the two configura¬ 
tions. 


4 Advantages of the vertical configuration 

Despite the larger number of UCNs that can be stored in the horizontal setup, 
there are distinct advantages for the vertical setup in terms of systematic effects. 

In the horizontal configuration, the two magnetic mirrors can trap the charged 
decay products. In-situ counting with a high efficiency is possible only by ex¬ 
tracting the protons with high-voltage electrodes (> 5kV) due to their low 
kinetic energy of up to 750 eV only. The trapped electrons can produce ions 
and electrons in the residual gas of the vacuum, and the presence of the electric 
field can accelerate them causing production of secondary charges up to electri¬ 
cal break down. For instance, even in the aSPECT neutron /3-decay apparatus, 
with its excellent vacuum conditions, this effect has been problematic[29, 30] . 
The vertical configuration avoids the electron trapping problem since the up¬ 
per coil is not required for UCN confinement. Moreover, this allows electron 
counting as the detection scheme and thus removes the need for the high-voltage 
system. 

A key systematic effect is poor cleaning of above-threshold UCNs from the 
trap. For a horizontal configuration, a scheme of ramping up and down the bias 
field has been used for this, however this results in a significant loss (~ 50%) 
of UCNs|3T|. The use of a UCN reflecting paddle along the length of a trap to 
induce mode-mixing reflections has also been suggested[35]. However, removing 
the paddle from the reach of UCNs after cleaning by rotation or by sliding it out 
of the trap causes undesired doppler heating of the UCNs (for the latter case, 
















heating can occur after non-specular reflections [33|). In the vertical, magneto- 
gravitational configuration, UCNs have to fall and reflect off the bottom of the 
trap. The insertion of a piston that induces mode-mixing reflections from the 
bottom can thus be used for cleaning. Furthermore, retracting it does not cause 
doppler heating. The details of this procedure will be published elsewhere. 

Finally, an advantage offered by this magnetic bottle is the large end open¬ 
ing area to volume ratio. When it is vertical, the opening at the bottom will 
allow the live monitoring of depolarized or warmed UCNs with a high efficiency. 
This will allow direct measurements of these key systematics. Furthermore, if 
a fill-and-empty scheme is used to determine r n —as will be the case of initial 
measurements—emptying of the trap from the bottom will have weak sensitiv¬ 
ity to the phase-space evolution of the UCNs, a key topic discussed in these 
proceedings. 


5 Conclusion 

A detailed description of the magnetic fields from the Halbach-octupole array 
and the superconducting coil assembly design has been given. The trapping 
potential and the effective volumes for a horizontal and a vertical geometry were 
compared, with the former allowing the storage of a factor ~ 2.5 more UCNs. 
However, the advantages for controlling key systematic effects offered by the 
vertical setup outweighs the decreased statistics, which is addressed with the 
development of the new, compact superfluid helium UCN source [25]. This has 
led to the magneto-gravitational design for our planned high-precision neutron 
lifetime measurements. 
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